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We study the resonant neutrinoless double-electron capture (Oi^ECEC) in ^^^Gd, ^"''Er and ^*'^W 
atoms, associated with the ground-state to ground-state nuclear transitions. The corresponding 
matrix elements are calculated within the deformed QRPA using the realistic Bonn-CD nucleon- 
nucleon interaction. The half-lives are estimated with the use of the most recent precision data 
on the Q-values of these processes. Perspectives of experimental search for the Or/ECEC with the 
isotopes Gd, ^'^''Er and ^*°W are discussed. 
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I. INTRODUCTION 

One of the unsolved mysteries of today's particle 
physics and cosmology is the question of whether neu- 
trinos are Dirac or Majorana particles. In the first case 
neutrinos and antineutrinos are fundamentally different, 
whereas in the second case neutrinos and antineutrinos 
are identical. Theoretical arguments in favor of Majo- 
rana neutrinos exist for decades in connection with the 
smallness of neutrino masses. Some grand unified the- 
ories explain smallness of the masses, e.g., in a seesaw 
scenario with heavy Majorana leptons 

Neutrinoless double-beta decay (Oj//3/3), 

{A,Z)^{A,Z + 2) + e- +e-, (1) 

is being considered as the unique practical tool for de- 
termining the nature of neutrinos (for a review, see 0). 
The most favorable decays for the experimental search 
are those with high mass difference between the ground 
state neutral atoms, i.e. Q-values, in which the parent 
nuclei decay to ground states of the daughter nuclei. 

The mere observation of neutrinoless double-electron 
capture (Oj/ECEC), 

e- +e- + iA,Z)^iA,Z-2r, (2) 

could also prove the Majorana nature of neutrinos as well 
as the violation of the total lepton number conservation. 
A double asterisk in Eq. ^ means that, in general, the 
final atom {A, Z — 2) is excited with respect to both the 
electron shell, due to formation of two vacancies for the 
electrons, and the nucleus. This process, as noted long 
time ago by Bernabeu, De Rujula and Jarlskog 0, may 



have a resonant character under the condition of degen- 
eracy of the masses of initial and intermediate atoms. 
In contrast to the Qi^/3f3 decays, in the Oi^ECEC capture 
small Q-values are favorable. The capture rate is a sen- 
sitive measure of the neutrino mass. 

The resonance enhancement can increase the probabil- 
ity of capture by many orders of magnitude. In searches 
for Majorana neutrinos, the neutrinoless double-electron 
capture can compete with the neutrinoless double-/3 de- 
cay provided the resonance condition is satisfied within 
a few tens of electron-volts. So far, however, there was 
no way to identify promising isotopes for experimental 
search of Oi^ECEC, because of poor experimental accu- 
racy of measurement of Q-values which until recently 
were known with uncertainties of 1 - 10 keV only [J]. 
Progress in precision measurement of atomic masses with 
Penning traps 0-0] has revived the interest in the old 
idea on the resonance Oz^ECEC capture. 

Sujkowski and Wycech 8] and Lukaszuk et al. @ ana- 
lyzed the Oi^ECEC process for nuclear 0+ O"*" tran- 
sitions accompanied by a photon emission in the res- 
onance and non-resonance modes. The physical back- 
ground for the process was calculated. A new theoreti- 
cal approach developed by Simkovic and Krivoruchenko 
and Krivoruchenko et al. [ll| allowed a unified de- 
scription of the oscillations of stable and quasistationary 
atoms, which take place with violation of the total lepton 
number conservation and are followed by de-excitation 
with emission of photons. Based on the most recent data 
and realistic evaluation of the decay half-lives, a com- 
plete list of the most perspective isotopes for which the 
Oj^ECEC capture ma y h ave the resonance enhancement 
was provided in Ref. |ll| for further experimental study. 
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Some isotopes such as -"^^^Dy have several closely-lying 
resonance levels. A more accurate measurement of Q- 
value of atoms ^^^Dy and ^^^Gd confirmed the existence 
of overlapping Oi'ECEC resonance levels [l^ . Assuming 
an effective mass for the Majorana neutrino of 50 meV 
and an appropriate value of nuclear matrix element, half- 
lives of some of the isotopes were found to be as low as 
10^^ years in the unitary limit, which is one order of 
magnitude shorter than the Oi^/3/3 half-life of ^®Ge for 
the same mass of Majorana neutrino [T3| . 

In high-Z atoms, the electrons in inner shells are mov- 
ing with relativistic velocities. Effects associated with 
the relativistic structure of the electron shells reduce the 
Oj^ECEC half-lives by almost one order of magnitude. 
In contrast to the non-relativistic theory, the capture of 
electrons from the npi/2 states is only moderately sup- 
pressed in comparison with the capture from the nsi/2 
states. In the relativistic formalism, selection rules ap- 
pear to require that nuclear transitions with a change in 
the nuclear spin A J > 2 are strongly suppressed. The rel- 
ativistic effects also enhance the violation of parity in the 
the Oi^ECEC process, as a result of which nuclear tran- 
sitions 0+ —> 0^,1^ become all attainable for a mixed 
capture of s- and p-wave electrons [ll| . A similar effect 
occurs due to weak right-hand currents as discussed by 
Vergados [l^ . 

Recently there has been fast progress in the measure- 
ment of atomic masses with the help of Penning traps. 
The accuracy of Q-values at around 100 eV was achieved 
[l^ [T5i-,23.] , which has already allowed to exclude a num- 
ber of isotopes from the list of the most promising candi- 
dates for searching the neutrinoless double-electron cap- 
ture. The best candidate is currently ^^^Gd, which al- 
though does not reach the unitary limit, however, un- 
dergoes a significant increase in the capture rate due to 
the proximity to the resonance level |20j . Further pre- 
cise measurements of masses of prospective isotopes are 
vigorously encouraged to continue. 

Neutrinoless double-electron capture has a number of 
important advantages with respect to experimental sig- 
natures and background conditions. In recent years, ex- 
perimental searches for the capture process were contin- 
ued [H-Hol. New upper limits of about 10^^ - 10^^ years 
for the half-lives of ''^Se, ^°^Cd and ^^^Sn have been ob- 
tained [23,[li,[2l[2l|30|. To make further progress new 
experimental data on excited states of finite-state nu- 
clei (excitation energy, angular momentum, parity) and 
precise calculations of transition nuclear matrix elements 
(NMEs) are required. Among the promising isotopes, 
^^^Gd, -'^^'^Er and ^^'^W have likely resonance transitions 
to the 0"'" ground states of the final nuclei. ^ 



1 The transition of ^°^Cd to an excited state of ^'^^Pd with the 
nuclear excitation energy of 2717.59 keV was examined in detail 
in Ref. l3l| for Jf = 0. In Ref. [Till it was noted that, as long 
as this level 7-decays by 100% into the 3"*" state at 1557.68 keV, 
the assignment Jj^ = is excluded. 



In this paper, accurate calculations of the Oi^ECEC 
half-lives of ^^^Gd, ^^^Et and ^^°W are performed. The 
electron wave functions in the atoms are treated in the 
relativistic Dirac-Hartree-Fock approximation [32*1 . The 
nuclear matrix elements are calculated within the proton- 
neutron deformed quasiparticle random-phase approxi- 
mation (deformed QRPA) with a realistic residual inter- 
action [33- 36|. 



II. RESONANT NEUTRINOLESS 
DOUBLE-ELECTRON CAPTURE 

The Oi^ECEC leads to a violation of conservation of 
total lepton number by two units. If the process is due to 
the Majorana neutrino exchange mechanism, the capture 
rate is determined by the effective Majorana neutrino 
mass 



TO/3/3 



(3) 



Here, Uei is the element of the Pontecorvo-Maki- 
Nakagawa-Sakata neutrino mixing matrix (stI . [38| and 
rrii are the diagonal Majorana neutrino masses. 
The half-hfe of the Oi/ECEC has the form 



nOi/ECEC 
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-^1/2 pOi^ECEC ' 

where the decay width is given by 
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with 



A = Ma,z - = Q-Ba. 



(4) 



(5) 



(6) 



being the difference of masses of the initial and final ex- 
cited atoms with masses Ma.z and 2,-11 respectively. 
Bab = Ea + Eb + Ec is the energy of two electron holes, 
whose quantum numbers (n, j, I) are denoted by indices a 
and b, Ec is the interaction energy of the two holes, while 
Tab is the width of the excited final atom with the elec- 
tron holes. Since we consider the transitions to ground 
states of the final nuclei, the single asterisk for 
refers to the excitation of electron shells only. 

Having factorized the electron shell structure and the 
nuclear matrix element, the lepton number violating 
transition amplitude can be represented as 11 1 



9\ 



Vab=m0pGj^{Fab)M' 



(7) 



Here, Gp ^ Gp cos 9c, where 9c is the Cabbibo angle, qa 
is the axial- vector coupling constant, R is the nuclear ra- 
dius, {Fab) is a combination of averaged upper and lower 
bispinor components of the atomic electron wave func- 
tions defined in Ref. [lH. M^" is the nuclear matrix el- 
ement of 0+ ground state to 0+ ground state transition, 
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TABLE I: Phenomenological paring gaps for protons Ap and neutrons A,i and deformation parameters /3 deduced from intrinsic 
quadrupole moments measured by the Coulomb excitation reorientation technique (/3qp , sign is given exphcitly if known) [41I ] 
and B{E2) values {Pb(E2)) [ill in ^^^Gd, ^''*Er and ^*"W. /32 is the deformation parameter of Woods-Saxon mean field fitted 
to reproduce the experimental quadrupole moments. {BCSi\BCS f) is the BCS overlap between the initial and final BCS 
vacua [4^. 

Initial (final) Ap A„ Pq^ Pb{E2) P2 {BCS^\BCSf) 

nucleus (MeV) (MeV) 

152q^ (i52gjn) 1.478 (1.117) 1.179 (1.192) (+0.29) 0.212 (0.306) 0.166 (0.256) 0.44 
i64j,j. (i64j-,y) (0.879) 1.020 (0.825) 0.36 (+0.32) 0.333 (0.348) 0.289 (0.302) 0.73 

180^ (isoHf) 0.927 (0.832) 0.788 (0.713) 0.27 (+0.27) 0.252 (0.273) 0.237 (0.244) 0.75 



which is a sum of the Fermi (F), Gamow- Teller (GT), 
and tensor (T) contributions [3^, : 



M°^ = -^+Mgt 



(8) 



In comparison with the corresponding 0^/3/3 decays, the 
isospin operators r+ of nucleons entering the NMEs are 
replaced by r^. 



III. CALCULATION OF NMES 



Nuclei participating in the Oz^ECEC ground state 



to ground state nuclear transitions ^"'^'^Gd ^^^Sm, 
i64ej. ^ I64j3y 180^ ^ isojjf defomied. The 

deformation parameter /3 = -y/Tr/SQj,/ {Zr'^) can be de- 
duced from the intrinsic quadrupole moment Qp of the 
first 2+ state measured by the Coulomb excitation re- 
orientation technique (rc is the root mean square charge 
radius). Unfortunately, the electric quadrupole moment 
of ^^^Gd has not been measured yet by this method [4l| . 
Alternatively, the deformation parameter /? can be ex- 
tracted from the values of measured E2 transition prob- 
ability, ( I Qp I = ^yl6^^B{E2)/5e'^, the sign cannot be ex- 
tracted [43 )• From Table H] one can see that the defor- 
mation parameters determined in both ways agree quite 
well with each other. 

NMEs of the considered Oi^ECEC transitions are calcu- 
lated within the deformed QRPA with a realistic nucleon- 
nucleon interaction [33l - [36| . The details of the formal- 
ism for the Oi//3/3-decay NME are given in Refs. [35l[36|. 
The needed generalization of the basic equations of 
Refs. [H m to the case of Oi^ECEC NME is straight- 
forwardly achieved by replacing r+ with t~ operator. 
We note that in the calculation of M'^^ within the de- 
formed QRPA the tensor contribution has been neglected 
till now. Within the spherical QRPA this component of 
M^" reduces its total value by less than 10 % [io, SI El. 

The single particle states are those of the axially sym- 
metric Woods-Saxon mean field and are expressed in the 
basis of an axially-deformed harmonic oscillator states. 
The parametrization of the mean field is adopted from 
the spherical calculations of Refs. [13, 113, EE]. The 
single-particle model space consists of 4 — 61iuj shells in 



the spherical limit. Only the quadrupole deformation is 
taken into account in the calculation. The fitted val- 
ues of the parameter /32 of the deformed Woods-Saxon 
mean field, which allow us to reproduce the experimen- 
tal I3b{E2), are shown in TablelH The spherical limit (i.e., 
/?2 = 0) is considered as well, to compare with the earlier 
results fl3| obtained in the spherical QRPA. 

We use the nuclear Brueckner G-matrix, obtained by a 
solution of the Bethe-Goldstone equation with the Bonn- 
CD one boson exchange nucleon-nucleon potential, as a 
residual interaction. The BCS equations are solved to 
obtain occupation amplitudes and quasiparticle energies, 
constituents of the nuclear Hamiltonian. The pairing in- 
teractions are adjusted to fit the empirical pairing gaps 
for protons and neutrons. The gap parameters are deter- 
mined phenomenologically from the odd-even mass dif- 
ferences through a symmetric five-term formula involving 
the experimental binding energies. The values obtained 
from this procedure for the nuclei under consideration 
can be seen in Table H] The calculated BCS overlap fac- 
tors of the initial and final quasiparticle mean fields are 
given in the last column of Table H] This factor repre- 
sents a possible suppression of the Oi'ECEC NME due to 
different deformations of initial and final nuclei [IJ . 

To solve the deformed QRPA equations, one has to fix 
the particle-hole gph and particle-particle gpp renormal- 
ization factors of the residual interaction of the nuclear 
Hamiltonian. Since experimental information about the 
position of the Gamow- Teller giant resonance for ^^^Gd, 
^^■^Er and ^^"W is not available we consider gph = 0.9 
like in previous calculation of the Oz^/3/3-decay NME of 
^^°Nd and ^^°Gd. In [13] it was proposed to adjust the 
particle-particle strength parameter gpp to the measured 
2:^/3/3-decay half-lives, i.e., to reproduce the experimen- 
tal value of the matrix element M^'ip. This procedure 
makes the 0;//?/3-decay NMEs essentially independent of 
the size of the single-particle basis and the nuclear struc- 
ture input. Due to a small Q-value, the half-life of 
the double-electron capture with emission of two neu- 
trinos (2Z/ECEC) of i52Gd, i64Er and i**°W is too large 
to be measurable. Therefore, the parameter gpp of the 
QRPA is fixed by the assumption that the matrix ele- 
ment of the 2^ECEC process lies within the range 
(0, 0.10) MeV-^ Recall that Mg'^ for double /3-decaying 
nuclei from the region i28,i30rpg^ ^^^Xe and ^^"^Nd does 
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not exceed the above range by assuming the weak-axial 
couphng constant to be unquenched (g^ = 1.269) or 
quenched {gA — 1-0). As it will be shown below this 
procedure of fixing gpp is not a significant source of un- 
certainty in the calculated Oi/ECEC half-lives. 



TABLE II: The matrix element M"" for the 2i/ECEC of 
^'^^Gd, ^^"Er and ^*°W calculated within the spherical and 
deformed QRPA with realistic nucleon-nucleon interaction 
(Bonn-CD potential) gl]. 



Nucleus 












[MeV-^] 


sph. 


def. 


def. 






QRPA 


QRPA (/32 


= 0) QRPA 


i^^Gd 


0.10 


7.59 


7.50 


3.23 




0.00 


7.21 




2.67 




0.10 


6.12 


7.20 


2.64 




0.00 


5.94 




2.27 




0.10 


5.79 


6.22 


2.05 




0.00 


5.56 




1.79 



In Table in the NME M^" for ^^^Gd, ^^-^Er and ^^Ow 
calculated within the spherical and deformed QRPA are 
presented. There is a qualitative agreement between the 
results of the spherical QRPA and the spherical limit of 
the deformed QRPA (/32 = 0).^ Since the deformed nuclei 
are described within the adiabatic Bohr-Mottclson ap- 
proximation, the spherical limit of the deformed QRPA 
should be taken with caution. The differences can be 
attributed to the fact that within the spherical QRPA 
[iol I45I I one-body and two-body matrix elements enter- 
ing the expressions for Mq";^ and M°^, respectively, are 
calculated with single-particle wave functions approxi- 
mated by the spherical harmonic oscillator ones, while 
realistic Woods-Saxon single-particle wave functions are 
used in the deformed QRPA ^i^. In addition, M^" 
obtained within the spherical QRPA contains also Mt 
contribution, which can reduce its value by up to 10%. 
The results in Table HI] indicate that the nuclear defor- 
mation decreases the value of M°'' by more than factor 
of 2-3. We note that the largest suppression of M^" due 
to deformation is realized for A — 180 nuclear system, 
where deformations of initial and final nuclei are compa- 
rable (same sign is assumed, see Table H]). It means that 
the suppression of M°'' can be associated with the large 
deformation of initial and final nuclei and large value of 
A. Before the effect of deformation on for nuclei with 
smaller A was associated with difference in deformations 
of the initial and final nuclei j46„ . 



IV. THE Oz/ECEC HALF-LIVES OF ^^^Gd, ^'^^Er 
AND i*°W 

Equations (4) - (7) imply that the half-lives are deter- 
mined by the following properties of the initial and final 
atoms: 

i) The mass difference A ^ determines the proximity 
to the resonance condition and ultimately the magnitude 
of the effect. It depends on the Q-value and the energy 
of two electron vacancies in the final atom. The selection 
rule of Ref. [ll| imply that it suffices to consider the 
capture of nsi^2 8'iid npi/2 electrons. 



TABLE III: Upper and lower average components of the 
Dirac bispinors (/) and {g) in ^^^Gdand ^''^Erfor lsi/2, 2si/2, 
3si/2, and 2pi/2 electron shells (in keV^'^^). The upper lines 
give solutions based on the Dirac equation in the Coulomb 
field [ll[ , the lower lines give the same solutions based on the 
Dirac-Hartree-Fock method [321 . 



Average 



'Gd 



= Er 



(/(IS1/2)) 1.33 X lO-* 
1.27 X W 



(/(2si/2)) 

(/(3Sl/2)) 

(ff(2pi/2)> 



5.20 X 10^ 
4.59 X 10^ 
2.84 X 10^ 
2.15 X 10^ 
-1.12 X 10^ 
-9.53 X 10^ 



1.57 X 10"* 
1.50 X 10" 
6.20 X 10^ 
5.46 X 10^ 
3.39 X 10^ 
2.60 X 10^ 
-1.43 X 10^ 
-1.22 X 10^ 



^ Note, thatOp/i = 0.9 is used here in all calculations while the 
results of [ijj were obtained with = 1.0. 



Recently, Q- values in the Oi'ECEC transitions of ^^^Gd 
[13], ^'^^Er 0, and i®°W [22] have been remeasured us- 
ing Penning-trap mass spectrometry with an uncertainty 
of 180 eV, 120 eV, and 270 eV, respectively. 

In atomic physics, the electron binding energies are 
usually measured to within a few eV. We used the bind- 
ing energies of single electron holes Ea from Ref. [isj . 
Noticeable corrections arise from the Coulomb interac- 
tion between two holes Ec- Coulomb interaction energy 
is calculated on the basis of the Dirac equation with ac- 
count of the screening of the nuclear charge, which gives 
an accuracy comparable to the present experimental er- 
rors in Q- values. More accurate estimates can be ob- 
tained by averaging the Fermi-Breit potential over the 
states of the two vacancies. 

ii) In the unitary limit the capture rate is inversely 
proportional to the width of the decay of excited atoms 
in the final state. The radiative width of the electron 
shell with two holes is estimated as Tab = Pa + Pfe on the 
basis of the measured and recommended values of radia- 
tive widths Pa of the single vacancies [i^ . The radiative 
width also determines the desired accuracy in measuring 
Q- values. 

iii) In nuclei with large values of Z, electrons in the 
lower shells are moving at a speed close to the speed 
of light. The wave functions of electrons in the Coulomb 
field plus the self-consistent field of neighboring electrons 
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TABLE IV: The calculated Oi^ECEC half- lives of ^^^^Gd, ^'^^Er, and ^^^''W for m^^j = 50 meV. The second and third columns 
show the quantum numbers of the electron holes. Here, n is the principal quantum number, j is the total angular momentum, 
and I is the orbital momentum. Shown in the columns four, five and six are the hole energies and their Coulomb interaction 
energy (in units of keV). The column seven shows the radiative widths of the excited electron shells. The column eight shows 
the mass difference of the initial and final atoms. The last two columns show the minimum and maximum half-lives of the 
O/^ECEC transitions (in years). The masses and the energies are in keV. 



Nucleus {n2jl)a {n2jl)t 


Ea 


Eb 


Ec 


Tab (kcV) 


A (kcV) 


Tr;2 (y) 


Tr/T (y) 


is^Gd 


110 


210 


46.83 


7.74 


0.34 


2.3 X 10~ 


-2 


-0.83 ±0.18 


4.7 X 10^* 


4.8 X 10^^ 




110 


211 


46.83 


7.31 


0.32 


2.3 X 10" 


-2 


-1.27 ±0.18 


4.2 X 10^^ 


1.1 X 10^^ 




110 


310 


46.83 


1.72 


0.11 


3.2 X 10" 


-2 


-7.07 ±0.18 


9.4 X 10^^ 


1.1 X 10^^ 


164 


210 


210 


9.05 


9.05 


0.22 


8.6 X 10" 


-3 


-6.82 ± 0.12 


7.5 X 10^^ 


8.4 X 10^^ 




210 


211 


9.05 


8.58 


0.23 


8.3 X 10" 


-3 


-7.28 ±0.12 


4.2 X 10^'' 


4.6 X 10^* 




210 


310 


9.05 


2.05 


0.11 


1.8 X 10" 


'2 


-13.92 ±0.12 


3.5 X 10^^ 


3.9 X 10^^ 




110 


110 


63.35 


63.35 


1.26 


7.2 X 10" 


-2 


-11.24 ± 0.27 


1.3 X 10^^ 


1.8 X 10^^ 



must be considered in a relativistic approach. The upper 
and lower radial functions of the Dirac bispinors are aver- 
aged over the volume of the nucleus. Relativistic effects 
lead to an increase in the probability of Oi^ECEC process. 
Estimates of the average values of the radial component 
of the Dirac wave functions inside the nucleus can be 
obtained using well-known wave functions of electrons 
in the Coulomb field, or by solving many-body problem 
using the Dirac- Hartree-Fock method. The nsi/2 elec- 
tron capture is proportional to the mean value of the 
upper bispinor radial component, while the capture from 
the npi/2 states is determined by the mean value of the 
lower radial components of the Dirac wave function. In 
Table Hill the average values of the upper and lower radial 
components of Ref. are compared with the values ob- 
tained in the Dirac- Hartree-Fock method [sl] for ^^^Gd 
and ^^^Er. The agreement is quite good. 

iv) The probability of capture is proportional to the 
square of the nuclear matrix element, which we discussed 
in the previous section in detail. 

In Table IIVI the maximum and minimum values of the 
Oi^ECEC half-lives (in years) of ^^^Gd, ^^^Ei and ^^°W 
for m^^ = 50 meV assuming the capture of most fa- 
vored atomic electrons, the associated values of the bind- 
ing energies and the Coulomb interaction energy of the 
two holes, and the mass difference A of initial and fi- 
nal excited atoms are shown. The estimates for T™™ 
are taken for the minimum mass difference deviating by 
not more than three standard errors from the experi- 
mental mean value unlike in [l3| , where the unitary limit 
was considered. Among the transitions, the favoured one 
is the capture of electrons from K and L shells in the 
case of ^^^Gd, which results in the half- life in the range 
4.7 X 10^^ ^4.8 X 10^'' years. This transition is still rather 
far from the resonant level. The half-life appears thereby 
2-3 orders of magnitude greater as compared to the 
halt life of 01^1313 decay of ^^Ge ;4^. 

V. CONCLUSION 

As shown in Ref. [Tl| , the Oi'ECEC half- lives of a dozen 
of isotopes are comparable to the shortest half-lives of 



the QvjSp decays of nuclei provided the resonance con- 
dition is matched with an accuracy of tens of electron- 
volts. Among the promising isotopes ^^^Gd, ^^''Er, and 
^*°W were found to be associated with the transitions 
between grotmd states of the nuclei. The estimates of 
the Oi/ECEC half-lives were recently improved by more 
accurate measurements of Q-values for ^^^Gd ^20,], ^^''Er 
p^ . and ^^"W 23] in a Penning trap. 

In this paper, we have made a further step to refine the 
estimates of the half-lives by going beyond the spheri- 
cal approximation in the calculation of nuclear matrix 
elements. We found within the deformed QRPA that 
the deformation of the nuclei leads to suppression of the 
NMEs by the factor of 2-3 as compared with the spher- 
ical limit. The suppression of NMEs depends not only 
on the relative deformation of the initial and final nuclei, 
but also on their absolute values. 

We conclude that the Oi^ECEC half-life of ^^^Gd is 2- 
3 orders of magnitude longer than the half-life of Qvfip 
decay of ^®Ge corresponding to the same value of the 
Majorana neutrino mass. Our calculation excludes ^^"'Er 
and ^®°W from the list of prospective isotopes to search 
for the neutrinoless double-electron capture. 
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